Patch coalescence as a mechanism for eukaryotic directional sensing 
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Eukaryotic cells possess a sensible chemical compass allowing them to orient toward sources of 
soluble chemicals. The extracellular chemical signal triggers separation of the cell membrane into 
two domains populated by different phospholipid molecules and oriented along the signal anisotropy. 
We propose a theory of this polarization process, which is articulated into subsequent stages of germ 
nucleation, patch coarsening and merging into a single domain. We find that the polarization time, 
te, depends on the anisotropy degree e through the power law cx e~^, and that in a cell of radius 
R there should exist a threshold value eth oc R~^ for the smallest detectable anisotropy. 
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The cells of multicellular organisms are endowed with 
a chemical compass of amazing sensitivity, formed as a 
result of billion years of evolution. Concentration differ- 
ences of the order of a few percent in the extracellular sol- 
uble attractant chemicals from side to side are sufficient 
to induce a chemical polarization of the membrane lead- 
ing to cell migration towards the signal source [ij. This 
way, a sensible amplifier of slight gradients in the distri- 
bution of chemicals in the surrounding environment is re- 
alized. Its relevance is easily understood if one recognizes 
that no multicellular organism could exist without the 
constituent cells being able of sensing directional signals. 
Directional sensing is actually essential both in embryo 
development, where tissue formation is realized through 
coordinated migration of specific cells guided by chemical 
signals, and in the adult organism, where chemical sig- 
nals guide white blood cells to the sites of inflammation 
and platelets to sites of wound repair. The main steps 
of the process are as follows (see the reviews Q and 0]). 
As a response to the attractant signal, the cell membrane 
is polarized, afterwards inducing differentiated polymer- 
ization of the cell cytoskeleton in its proximity. The re- 
sulting unbalance, triggered by a well characterized cas- 
cade of chemical reactions, leads to the formation of a 
growing head and a retracting tail, in such a way that 
the cell starts to drift towards the source of the signal. 
The initial part of this process is constituted by the early 
chemical polarization of the cell membrane. In this letter 
we propose a simple phenomenological scheme providing 
a universal description of this fundamental phenomenon. 

Membrane polarization can be recognized as a self- 
organization process governed by a network of diffusion- 
controlled chemical reactions. It is known that reaction- 
diffusion networks may become bistable in the presence 
of chemical feedback loops 0, ■ In spatially extended 
systems bistability may lead to the formation of compet- 
ing phases and to a phenomenology typical of first order 
phase transitions, such as metastability, nucleation and 
coarsening [1,0]. The polarized membrane state observed 
during directional sensing can therefore be interpreted as 



the coexistence of domains of two different phases. 

Let us briefly describe the chemical reactions which are 
responsible for directional sensing. The chemical factors 
clustering in complementary membrane domains are the 
phospholipids PIP2 and PIP3. Two enzymes, PI3K and 
PTEN, respectively transform PIP2 into PIP3 and vice- 
versa. The phospholipids are permanently bound to the 
inner face of the cell membrane, while PI3K and PTEN 
diffuse in the cell volume and are active only when they 
are adsorbed by the membrane. PI3K adsorption takes 
place through binding to receptors activated by the ex- 
tracellular attractant signal. This way, the external at- 
tractant field is coupled to the inner dynamic of the cell. 
PTEN adsorption takes place through binding to the 
PTEN product, PIP2. This process introduces a posi- 
tive feedback loop in the system dynamics [1, 0] . When 
the cell is not stimulated by an attractant signal the cell 
membrane is uniformly populated by PTEN and PIP2 
molecules. When a uniform receptor stimulation of a 
suitable amplitude is switched on, PI3K molecules bind 
to the membrane and shift its chemical balance toward a 
PIPa-rich phase, while PTEN desorbs. PIPa-rich germs 
are then nucleated in the PIP2-rich sea and PIPs-rich 
regions start to coexist with PIP2-rich ones (lo| . 

Two different regimes of cell polarization may be dis- 
tinguished. Anisotropy driven polarization induced by 
the presence of an attractant gradient is realized in a time 
of the order of a few minutes, and results in the forma- 
tion of a PIPa-rich domain on the membrane side closer 
to the attractant source and of a PIP2-rich domain in the 
complementary region [sl, Q. On the other hand, cells ex- 
posed to uniform distributions of an attractant polarize 
in random directions, in times of the order of an hour 
(see 6.17. [ll|). The existence of two clearly separated po- 
larization regimes is confirmed by the recent observation 
of a sensitivity threshold of the order of a few percent 
difference in the attractant molecule concentration from 
side to side [ij. Direct observation of the polarization 
process [1, Fig. 7a] implies the bound < 5 s for PTEN 
diffusion time in the cell volume, which is therefore much 
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less than the polarization time for both regimes. In this 
process, the amplitude of the cell stimulation is of cru- 
cial importance. At very low stimulation levels PTEN is 
not desorbed in a significant amount and no directional 
sensing takes place. At very high stimulation levels a 
homogeneous PIPs-rich phase is realized and directional 
sensing again does not take place. There exists therefore 
an optimal attractant concentration, such that below it 
the minority phase is PIPs-rich and above it is PIP2-rich. 

Numerical simulations of the directional sensing net- 
work performed with the use of realistic physical and 
kinetic parameters have shown that under appropriate 
conditions the biochemical network is indeed bistable, 
and that it undergoes spontaneous separation in chem- 
ically different phases, rich in PIP2 and PIP3, respec- 
tively 0, [13]. A 5% anisotropic component in the cell 
stimulation accelerates cell polarization and correspond- 
ingly decreases the characteristic time needed for com- 
plete phase separation by more than one order of mag- 
nitude: fast, anisotropy driven polarization is realized in 
times of the order of a minute, while slow, stochastic po- 
larization is realized in times of the order of one hour, in 
accordance with experimentally observed times. In the 
numerical experiments, when PIP2 is the minority phase, 
the evolution leading to phase separation consists of an 
early nucleation regime, resulting in the formation of iso- 
lated PIP2-rich patches and a late coarsening process, 
where large patches of the PIP2-rich phase grow at the 
expense of the evaporation of smaller ones, similarly to 
what happens in the case of first order phase transitions 
in a liquid-gas system or in the precipitation of a su- 
persaturated solution • Finally, the patches condense 
into a single large cluster, leading to a stationary state 
characterized by the coexistence of a PIP2- and a PIP3- 
rich domain. However, the dynamics of the directional 
sensing network differs from that of otherwise similar pro- 
cesses, such as the precipitation of a supersaturated so- 
lution, under one important respect. When precipita- 
tion nuclei in a supersaturated solution dissolve, matter 
is transferred to larger nuclei through diffusion in the 
surrounding medium. In contrast, in the directional sens- 
ing network enzyme-substrate patches evaporate through 
desorption of the PTEN enzyme from the membrane, 
which is then transferred to other patches through dif- 
fusion in the cell volume. Therefore the transformation 
of PIP3 into PIP2 molecules cannot be described at the 
membrane level as a local, diffusion-like process as is the 
case with the adsorption-desorption process from precip- 
itation nuclei; in particular, there is no local conservation 
of the number of PIP2 molecules. 

The above summarized scenario can be put on a firm 
analytical ground resorting to the kinetic theory of first 
order phase transitions [71, llJ, llJ]. In this theory, af- 
ter germ nucleation larger patches of the stable phase 
grow at the expense of smaller patches which shrink, 
leading to scaling laws and universal probability distri- 



bution of patch sizes. We shall now show how the ideas 
of the Lifshitz-Slyozov theory may be adapted to 
our problem to deduce simple scaling laws for the mem- 
brane polarization time and explain most of the observed 
phenomenology. We discuss here the case when PIP 2 is 
the minority phase (the other case being symmetric). In 
this case PIP2-rich patches are formed inside the PIP3- 
rich sea (see Fig. [1^). We restrict our consideration to 
approximately circular patches of the PIP2-rich phase, 
which are expected to dominate over different geome- 
tries due to the presence of a linear tension between the 
two phases. The free energy of a PIP2-rich patch of ra- 
dius a can be written on phenomenological grounds as 
F = — TT "00^ -|- 27r cr a, where a is the linear tension of 
the interface with the surrounding PIPs-rich phase and 
■0 represents the degree of metastability 13|, which is a 
function of the concentration of PTEN molecules in the 
cell volume and of the concentration of extracellular at- 
tractant. 

According to the kinetic theory of first order phase 
transitions, the equation of growth of a patch is dissi- 
pative. In the absence of a local conservation law the 
equation for a circular patch can be written as F dta = 
—dF/da, where F(a) is a damping coefficient Q. Since 
energy dissipation occurs mainly along the perimeter of 
the interface between the two phases, F may be written 
as 27ra7, where 7 is a constant, and we get 



^dta = -ip - a /a + ^ , 



(1) 



where the noise term ^ represents thermal fluctuations. 
The fluctuations are responsible for the formation of an 
initial population of patches with varying radii a [3, §99] . 
Patches with a smaller than the critical radius Qc = cr/tp 
are mainly dissolved while most patches with a > Oc sur- 
vive and grow because of the gain in free energy. At 
initial time, Oc is of the order of the thickness oq of 
the interface between the two phases 0, [l^. As long 
as the area occupied by patches of the PIP2-rich phase 
grows, the degree of metastability ip decreases, some of 
the patches that were initially growing become undercrit- 
ical and shrink, large patches start "feeding" on smaller 
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ones, and the total number of patches diminishes 
§100]. In the final stage of this process a single domain of 
the PIP2-rich phase is formed coexisting with the PIP3- 
rich phase, see Fig. [!};. However, the details of the pro- 
cess leading to this final stage depend on the external 
conditions, and, particularly, on the degree of anisotropy 
of the attractant signal. 

The population of patches can be described in terms 
of the size distribution function /(a) such that /(a)Aa 
gives the number of patches with sizes in the interval 
(a, a + Aa). An important simplification comes from the 
fact that for patches with a > Oc the noise term ^ in ([1]) 
becomes negligible [l^l. This means that the stochas- 
tic nature of the problem enters mainly in the formation 
of the initial distribution of patch sizes /(a), while for 
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a) t <te b) t ~ te 



2) t>U 



FIG. 1: Patch growth in the presence of a shght gradient of 
attractant activation directed from left to right. The PIP3- 
and PIP2-rich phases are respectively light and dark grey. 



a > flc the time evolution of /(a) is dictated by the de- 
terministic part of ([1]), from which the kinetic equation 
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follows 0, [l^- With the chosen normalization J /(a) da 
represents the total number of patches, a quantity which 
is monotonically decreasing in time according to the pre- 
viously described "coarsening" dynamics [14|. Eg. (|2l) is 
valid as long as / /(a) da is much larger than 1 |13j . 

To obtain a closed system of equations we need an ad- 
ditional equation for the time evolution of the metasta- 
bility degree ij) . In the case of isotropic stimulation ij: 
does not depend on the position on the membrane and is 
instead only a function of time. Since diffusion of PTEN 
molecules in the cell volume is faster than phospholipid 
diffusion on the membrane we can regard the concen- 
tration of PTEN molecules in the volume as uniform ^ . 
Moreover, fast PTEN diffusion also implies that -0 instan- 
taneously adjusts to the changes in the size distribution 
function. While the total number of patches diminishes 
as an effect of the coarsening dynamics, the total area 
occupied by the patches, as well as the total number of 
PIP2 molecules found in the patches, monotonically in- 
creases towards their respective equilibrium values. The 
metastability degree if) is equal to zero in equilibrium, 
and tends to zero in accordance with 



■0 oc j4 — / da TT f{t,a), 



(3) 



as the total patch area tends to its limit value A. Eq. ([3|) 
reflects the fact that in the asymptotic region ip is pro- 
portional to the excess concentration of PTEN molecules 
in the volume with respect to the equilibrium value, and 
therefore to the difference between the area occupied by 
the PIP2-rich phase at equilibrium and at current time. 
The law ([3]) is valid for t ^ to, where tg is the character- 
istic time needed for the formation of a germ of an alter- 
native phase that can be estimated asto '-^ ^ 0,0 0jll3|- 
Asymptotically, (|2l3p lead to the self-similar solution 



(4) 



where g{^) = ^(1 - (.)~^ exp[-2{l - ^)~^] if ^ < 1 
and g{^) — if ^ > 1. Similarly to what happens in 



Lifshitz-Slyozov theory 1J| the total number of patches 
decreases in time due to the evaporation of small patches: 
/ da /(a) oc t^^, and from (j4|) one gets (a) = a/ip = ac- 
The evolution of the size distribution function / governed 
by (dl stops at times of order i^,, where t^, is defined as 
the instant when the average patch size (a) reaches the 
cell size R. From the scaling law (a) oc we 
get ~ (i?/ao)^io- Eventually, at i ~ i*, a single 
PIP2-rich patch survives. Its orientation is determined 
by the random unbalance in the initial germ distribu- 
tion. Notice that in this derivation, following the lines 
of 0, m, 14 1, isotropy was essential to assume that ^ 
was uniform along the whole membrane surface. 

Let us now consider the case of an inhomogeneous acti- 
vation pattern. The inhomogeneity of the concentration 
distribution modifies the degree of metastability, which 
becomes a function of the position on the membrane sur- 
face. Since the distribution of PTEN molecules in the cell 
volume is homogeneous, it influences only the isotropic 
part of the metastability degree tp, which is a function of 
time, as previously. In contrast, the anisotropic part of 
the metastability degree, Sip, related to the external at- 
tractant inhomogeneity, does not depend on time. If the 
cell membrane has a nearly spherical form and a radius 
R much smaller than the characteristic scale of the ex- 
tracellular attractant distribution, then 6ip = — e ipo cos 9. 
Here ijjQ = (j/ao is the initial metastability degree, e is 
a dimensionless factor measuring the initial anisotropy 
degree, and 9 is the azimuthal angle on the cell surface. 
This way we obtain the equation 



7 dta = ip — etpo cos 9 — a/a + ^, 



(5) 



generahzing H]). As long as etpo <C ip, the first stage of 
patch growth proceeds approximately as in the isotropic 
case and ip decreases as t~^/^. However, at a time of or- 
der te, where is defined by the equation tp{t^) = eipo^ 
the perturbation e ipo cos 9 becomes comparable to ip and 
the process of polarization becomes anisotropic, so that 
patches in different regions get different average sizes, 
see Fig. [T}). From the scaling law ^ for tp one gets 

~ toe~^. For t > t^, the leading term in ([5]) becomes 
the perturbation eipocos9, implying that in the region 
closer to the source of the stimulation (cos 9 > 0) the 
PIP2-rich phase evaporates in a time which is easily esti- 
mated as being again of order , leading to the formation 
of a single PIP2-rich patch in the region further from the 
source of the stimulation (cos 9 < 0) and realizing com- 
plete polarization, as shown in Fig. [1^. 

The above scheme is valid as soon as the initial nucle- 
ation time to is significantly smaller than t^, an assump- 
tion which is compatible with the results of numerical 
experiments [j^]. On the other hand, the second stage of 
patch evolution occurs only if ^ t^- Otherwise, the 
presence of a gradient of attractant becomes irrelevant 
and only the stage of isotropic patch growth actually oc- 
curs. This condition implies that a smallest detectable 
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gradient exists, such that directional sensing is impossible 
below it. The threshold value eth for e is found by letting 

= ig. Since the product ipUc is a time-independent 
constant, we can simply compare its value at initial and 
final time when e = eth, obtaining eth = clq/R, which 
gives us the expression for the threshold anisotropy. 

It is interesting to estimate ao, and, consequently, eth, 
in terms of observable parameters. Comparing the char- 
acteristic patch surface and perimeter energy as a func- 
tion of the phospholipid diffusion coefficient D, surface 
phospholipid concentration c, surface concentration of ac- 
tivated receptors h, and the characteristic catalytic time 
T, one gets oq ~ {D t c/hY/"^ . Using parameter values 
from Ref. one gets oq ~ 1 /xm and eth ~ 10%. The 
value for etv, is compatible with the observations (the data 
from Ref. 'l] imply eth — 7% for dictyostelium) . 

One may wonder whether a cell may become polarized 
by the anisotropy produced by a spontaneous fluctuation 
in the extracellular distribution of attractant molecules 
or fluctuations in receptor- ligand binding Since eu- 
karyotic cells typically carry 10'*-10''^ receptors for at- 
tractant factors, one expects spontaneous fluctuations in 
the fraction of activated receptors to be of the order of 
10^^, a value which is comparable to observed anisotropy 
thresholds. However, to actually produce directed polar- 
ization the fluctuation should sustain itself for several 
minutes, i.e. for a time comparable to the characteris- 
tic polarization time. Such an event has very low prob- 
ability of being observed since the correlation time of 
the fluctuations determined by attractant diffusion at the 
cell scale and the characteristic times of receptor-ligand 
kinetics are much less than the polarization time. In- 
deed, the diffusion time is '--^ 1 s at the typical cell size 
10 /im, and the characteristic times of receptor-ligand ki- 
netics are also ~ 1 s (see online supporting information to 
Ref. [1]). Therefore, the direction of cell polarization in 
the case of a homogeneous distribution of attractant can 
only be determined by the inhomogeneity in the initial 
distribution of the positions of PIP2-rich germs produced 
by thermal fiuctuations. 

In conclusion, we have constructed a universal phe- 
nomenological description of the mechanism of direc- 
tional sensing in the eukaryotes based on the process of 
patch coarsening. This description implies the existence 
of two clearly separated polarization regimes depending 
on the presence or absence of an anisotropic component 
in the activation pattern produced by the extracellu- 
lar attractant factor, and the existence of a sensitivity 
threshold for the anisotropic component. Both results 
are in reasonable agreement with experimental observa- 
tions. Moreover, we predict that directed polarization 
time should scale as the inverse square of the relative 
signal anisotropy, a law that should be verifiable by di- 
rect observation. Our picture suggests that directed and 
stochastic polarization share a common mechanism, and 
that stochastic polarization should be the result of noise 



in subcellular and not in extracellular dynamics. Im- 
portantly, our picture does not depend on the details of 
the reactions involved, but only on the general structure 
of the directional sensing network and on its bistabil- 
ity. This means that the picture is robust not only with 
respect to variations of the kinetic and physical parame- 
ters, but also with respect to the identity of the chemical 
species involved. Indeed, PI3K and PTEN could be sub- 
stituted by, or synergize with, molecules endowed with 
similar enzymatic activity. An interesting speculation 
is that the bound eth = ao/R may explain why spatial 
directional sensing was developed only in the large eu- 
karyotic cells and not in smaller prokaryotes, whose di- 
rectional sensing mechanisms rely instead on the mea- 
surement of temporal variations in concentration gradi- 



ents 15| . Our bound derives from the intrinsic properties 
of polarization dynamics and is independent of the size 
criterion formulated in Ref. [l3|- The experimental ob- 
servation of selforganizcd phospholipid patches [I^ fol- 
lowing uniform attractant stimulation provides an initial 
confirmation of the validity of our scheme. To check the 
predictions of our theory, similar observations should be 
performed for the longer times characteristic of random 
and directed polarization, both under uniform attrac- 
tant activation and in the presence of accurately con- 
trolled concentration gradients. Experimental modula- 
tion of PTEN levels could be used to modify the overall 
size of patches and eventually switch off the patch for- 
mation mechanism. 
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